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ABSTRACT 



We propose a classical constrained Hamiltonian theory for 
the spin. After the Dirac treatment we show that due to the 
existence of second class constraints the Dirac brackets of the 
proposed theory represent the commutation relations for the spin. 
We show that the corresponding partition function, obtained via 
the Fadeev-Senjanovic procedure, coincides with the one obtained 
using coherent states. We also evaluate this partition function for 
the case of a single spin in a magnetic field. 
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1 Introduction 



The Heisenberg Model has been used for many years as a canonical model 
of magnetism in solids. After the discovery of High Tc superconductivity a 
renewed interest in this model emerged. This is principally due to the con- 
nection between the antiferromagnetic Heisenberg model and other models 
(Hubbard and t — J models) |1J which are important to describe the electronic 
structure of the High Tc superconductors. 

Although the Heisenberg model is quadratic in the spin operators it is 
highly non-trivial. What turns difficult the solution of the model are the 
complicated commutation relations among the spin components. This kind 
of commutation relations are indeed characteristic of systems described by 
Hubbard operators ||. In particular the spin algebra (su(2)) is the bosonic 
subalgebra of that spanned by the Hubbard operators in the so-called atomic 
representation of the t — J model, Osp(2, 2) [§]. 

When one is faced to a system with this commutation rules one usually 
implements some kind of decoupling using slave particles [[| . This new repre- 
sentation is formally exact and solves the above mentioned difficulty because 
the basic variables are conventional bosons and fermions. In particular a 
path-integral treatment is standard. However, in general, no exact calcula- 
tion can be carried out thoroughly and some kind of approximation is needed 
(V.gr. mean field, large s, large N, etc.). A natural question which emerges 
in this framework (first decoupling and then approximating) is what is an 
artifact of the decoupling and what is not (see for example || for a recent 
discussion). 

For the analysis of this question, a treatment of spin systems without the 
introduction of these fictitious slave particles would be of some importance. 
In particular, a path integral representation is a suitable formulation for this 
analysis. Indeed, such a path integral representation can be written using 
the coherent state method (also known as geometric quantization) || |(| . 

Although path integral methods have shown to be very powerful in vari- 
ous areas, the coherent state method is not a familiar subject in solid state 
physics. Moreover, as pointed in the literature 0,|6|, [7|] , there are some math- 
ematical subtleties in the derivation of this path integral representation that 
prevent for a complete confidence on the results. 

In this paper we present a constrained Hamiltonian system which after 
the Dirac treatment || reproduce the physics of the spin. Then we write 
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down its Feynman path integral following the Fadeev-Senjanovic method 
for constrained Hamiltonian systems ||. Our derivation is formally exact, 
independent of the coherent state method and valid irrespectively of the value 
of the spin s. Moreover, it has the virtue of setting the study of spin systems 
on canonical grounds. This is an important point in view of possible further 
applications of the formalism for spin 1/2 systems. 

Finally we explicitly calculate the free energy of a particle with fixed spin 
in an external magnetic field. The result gives an independent check for the 
validity of the expression of the partition function for all values of the spin. 
Besides, we think that the methods we use are very instructive and could be 
useful to develope approximation schemes for the Heisenberg model. 

The paper is organized as follows: in section II we present the constrained 
Hamiltonian system and we perform the Dirac treatment. In section III we 
find the path integral representation for the quantum mechanics and the 
partition function. We calculate the partition function of one spin in Section 
IV. Finally a brief summary of our main results is given in section V. 

2 The constrained Hamiltonian system 

The starting point for our analysis is the following Hamiltonian 

H = H(n) (1) 
and the additional set of primary constraints 

T x =P 1 -sA 1 {tS)kQ, (2) 

r 2 = P 2 - sA 2 (n) « 0, (3) 

r 3 = P 3 - sA 3 (n) « 0, (4) 
where P is the canonical momentum associated to the vector n and 

r 4 =|n| 2 -i^o (5) 
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Equations (1-5) define a classical constrained Hamiltonian system in the 
Dirac sense ||. Note that the constraints (2-4) are imposed in the weak 
sense. 

The classical Hamiltonian is a general function of a vector n. To simplify 
the notation we do not include lattice indices in the vector n but the following 
analysis is valid in the case where the spins live over a lattice. The Hamil- 
tonian (1) is independent of the momentum P associated to the vector field 
n. The corresponding momentum appears through the constraints (2-4). In 
these constraints we also have a general vector function A(n). In (2-4) s is 
an arbitrary number. It will become the size of the spin after quantization. 
Finally, the constraint (5) indicate the normalization of the vector n. 

According to the Dirac procedure we have to define the following total 
Hamiltonian: 

H T = H(n)+J2^i (6) 

i 

where the Lagrangian multipliers Aj are general functions of the coordinates 
and momenta. 

The consistency conditions of the theory require the preservation in time 
of the constraints (2-5) 

r i = {r i ,fz r } PB «o (7) 

In (7) {, }pb means the conventional Poisson Brackets. It is straightfor- 
ward to show that the system (7) determines uniquely the Lagrange multipli- 
ers (the explicit form of these multipliers is not necessary for the rest of the 
paper). Therefore the constraints are preserved in time and secondary con- 
straints are not generated. Moreover our theory is completely determined by 
the first class Hamiltonian (6) and the second class constraints (2-5) (there 
is no hidden gauge invariance in our system). 

Due to the second class character of the constraints, the next step is the 
calculation of the Dirac brackets (DB). For any two classical quantities A, 
B, the Dirac bracket is defined by 

{A, B} DB = {A, B} PB - {A, T^PBA^iTj, B} PB , (8) 
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where 



/\i j {Y j ,Y k }p B = 5 ik . (9) 

Using the PB, it can be shown that the antisymmetric matrix A" 1 has 
elements 



A 
A 
A 



-1 

'12 


JA 2 


dA 1 


dni 


dn 2 


-1 
13 


s\ dAs 
dni 


dA 1 
dn 3 


-1 

■14 


= -2m 




-1 
23 


,dA 3 
= S[ dn 2 


dA 2 
dn 3 


-1 
24 


= -2n 2 




-1 
34 


= ~2n 3 





A 2 

A,^ 1 = -2n 3 (10) 

In summary, we have a theory which is defined by a first class Hamiltonian 
(6) without arbitrary coefficients and the set of second class constraints (2)- 
(5). 

At this point the constraints are imposed as strong equations between 
coordinates and momenta. Also the PB are replaced by the DB. 



3 Quantization and Statistical Mechanics 

The canonical quantization can be now performed in the standard way. The 
Hamiltonian (6) have to be considered as the Hamiltonian operator. The 
constraints (2-5) become strong equations between operators and the com- 
mutation relations are obtained from the DB. The commutation relations 
between two operators 0\ and 2 are given by: 

[6i,6 2 ] = th{o 1 ,o 2 } DB . (n) 

The vector field A(n) is not specified yet. If we choose it as: 
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A A(n) = ft, (12) 
the DB between the components of ft reads: 



{ni,nj} DB = -j-. (13) 

In the last expression ijk are the integers 123 and its cyclic permutations. 
If we now identify shi — > Si, we obtain: 



[S l ,S J ] = ihS k (14) 

which corresponds to the commutation relations of the components of the 
spin operator. Then, we have proposed a classical constrained theory which 
after the Dirac procedure reproduce the commutation relations for the spin. 

Now, we are able to quantize the system using Feyman path integral. We 
shall use the method developed for constrained systems in |J. In this method 
the probability amplitude of the system which was at no at t — and will be 
at Tlx at time T can be written as: 



< n | n t >= f dftdP(^{8{Ti)){detI\- l ) 1/2 exp(^ C dt[Pft - H T \) (15) 

J ■ lb J 

Using the fact that det(A) is a constant independent of the fields (see 
(10)) and integrating over P, we obtain 

< n \ ni >= / dn5{\ ft | 2 — l)exp(— / dt[sA(n)n — H(ft)]). (16) 

J h Jo 

To obtain the partition function, we must carry out the integration over 
all periodic paths, make the change it — r and identify iT/h with the inverse 
of the temperature (3 Jl0| . Thus we have 



r 1 r0f>- -> df) 

Z= dn5{\ n | 2 -l)exp(-- dr[-sA{n)— + H{ft)\). (17) 

J periodic TlJO UT 
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Equations (16) and (17), with A given by (12), are exactly the expressions 
for the path integral for the spin system obtained via coherent states 0. 

Before going to the next section it is useful to make a few remarks about 
the reason for choosing the set of constrains (2)- (5). Our principal aim was to 
find a classical system containing the nature of the spin. To get a chance to 
find that, the classical system must contain constraints, and these constraints 
must be second class. The set of constraints (2)- (5) are arbitrary functions 
of the coordinates and momenta and our previous results showed that af- 
ter imposing condition (12) the quantized theory corresponds to the spin 
system. On the other hand, the classical constrained system we have intro- 
duced describe a system of massless particles living on unitary spheres. Each 
particle feels a magnetic field represented by the vector potential A. This 
vector potential corresponds to the one of a monopole placed at the center 
of the sphere as is stressed by condition (12). This mechanical analogy was 
previously suggested from the path integral expression obtained via coherent 
states || and has indeed inspired the election of our set of constraints. It is 
important to point out that once the constrained Hamiltonian system was 
defined we were able to show that the path integral representation is valid 
for any parameter s. This is a point not airways clear in the literature. 

4 An example 

Our starting point is eq. (|T7j) , with the Hamiltonian 

H{n) = -sn.B. (18) 

The first term in the action of eq.(|TTD can be written, with the aid of eq.(|T2]) 
and Stoke's theorem, as 

[P _ f f _ _ fP f 1 _ _ _ 

I dtA{n).d t n= \ \ n(t,r).da = dt drn.d t n x n. (19) 
Jo J J Jo Jo 

n(t, t) is a parametrization of the surface on S 2 enclosed by the trajectory 
n(t), that we chose so that n(t, 0) = n(t) and n(t, 1) = k. We then see 
that this term has a geometrical interpretation, which corresponds to the 
oriented area on the sphere S 2 enclosed by n(t). Terms of this kind are 
known as Wess-Zumino terms. 
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Using the explicit parametrization 
n(t,r) = 

(sin[(l - r)9(t)} cos <f>(t), sin[(l - r)9(t)} sin0(t), cos[(l - r)0(*)])(2O) 

where t 6 [0, r G [0, 1], the surface element is then given by 

n.d t n xn = <p(t)9(t) sin[(l - r)9(t)}. (21) 

The use of spherical coordinates for n allows for a straightforward imple- 
mentation of the constraint n 2 = 1 and the functional measure is 

Vn(t)5(n 2 - 1) = T>if(t)T>(cos9(t)). (22) 

Making use of the coordinates fl20|) , the partition function (17) can be written 

as 

Z = Vif(t)Vu(t)exp (is dt<p(t) (1 - u(t)) - sB u(t)dt\ , (23) 

where u(t) = cos(9(t)), the integral over r has been carried out in the WZ 
term and the functional integration is restricted to continuous closed trajec- 
tories. The corresponding boundary conditions are 

<p(P) = <p(0) + 2fcvr, 9(13) = 9(0) + 2k'it (24) 

and 

f((3) =f(0) + (2k + l)n, 9((3) = -9(0) + 2k'n. (25) 

We will first integrate out <f(t). It is convenient to integrate by parts the 
first term in the exponent of (^). Then it reads 

rP 

is (ip(p) - if(0)) (1 - u(0)) + is / (pit) u(t)dt. (26) 

J o 

The change of if along the trajectory, (A<f = <f(/3) — f(0)), is a constant 2kn 
and (2k + 1)tt for trajectories satisfying (PH) and respectively. 

The functional integral over f is now straightforward and the result is 
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J Vip{t)exp (is dt ip{t) (1 - u(t))j = ^ S[u(t)]exp (isAcp(l - u )) . 

(27) 

The ^-functional restricts the integral over u(t) to constant trajectories. 
For this reason any trajectory satisfying the boundary conditions fl25|) is 
wiped out from the partition function. The functional integral over the re- 
maining trajectories can be written as a sum of ordinary integrals, one for 
each posible value of Aip. 

Z = / du^exp (is2kn(l — u ) — sB(3uq) = 
k J ~ l 

J du J2^(s(l-u ) -n)exp(-sBp ). (28) 



The second equality follows from the Fourier representation of the 5-function. 
Integrating now over uq gives 

2s 

Z = Y,exp(-pB(s -r)), (29) 


which can be rewritten as 

s 

Z= exp(-pBm). (30) 

m=—s 

This expression for the partition function coincides with the well-known 
result obtained in elementary quantum statistical mechanics. 

It should be stressed that our result was obtained without the use of 
any approximation within the path-integral approach. One sould recall at 
this point that one of the mathematical subtleties involved in the formal 
construction of path-integrals is the restriction to continuous trajectories. 
Our computation, which explicitly uses this restriction, provides a check 
that this procedure leads to the correct result. 

It is important to point out that at the classical level there is no reason 
to choose the parameter s to be integer or half-integer. In this section we 
have found that this condition emerges naturally when the path integral for 
a 1-site problem is solved. 
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5 Summary 



In this paper we have presented a constrainted Hamiltonian theory for spin 
systems. Using the Dirac theory we have shown that our classical system has 
DB which can be identified with the commutation relation for the components 
of the spin. Following the recipe for quantizing constrained theories we have 
obtained a path integral representation for the spin system, which coincides 
with the expression obtained using coherent states 0. 

Our approach shows in a clear way that the path integral for the spin 
systems is free of any approximation. 

Finally, we have also recovered the partition function for a single spin in 
a magnetic field from our path integral expression. 

Natural extensions of the present work would be to study the case of 
Osp(2, 2) (related to the t-J model) and SU (3) groups. In fact, we are work- 
ing actually on this problem, but so far we have found that a generalization 
of this work is not straightforward and deserves a separate treatment. 
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